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The information about the effects of induced earth currents that can be obtained from analyses 
of observations and from calculations based on electromagnetic induction theory is examined. The 
nature of these effects is discussed and estimates of their magnitude for fields of varying extent and 
frequency are obtained. The influence of underlying geology and the effects of the oceans are con- 
sidered. The type of mathematical problem that needs to be solved to get more detailed information 
is described. 



1. Introduction 

It is well known that geomagnetic micropulsations 
of all types arise from sources external to the earth's 
surface. The primary oscillating electromagnetic 
fields of these external sources give rise, however, to 
induced electric currents within the earth, which con- 
siderably modify the observed surface fields. In 
seeking to discover the precise nature of the external 
sources, one must try to assess the effects of the in- 
duced earth currents and remove them from the 
surface observations. In this paper we consider 
what is known about these effects, and what can be 
inferred from analyses of the observations together 
with calculations based on electromagnetic induction 
theory. It is well known that the induced currents 
tend to reduce the amplitude of the Z-component and 
increase that of the //-component. An accurate 
estimate of these changes is needed, for example, 
when estimating the height and/or horizontal position 
of certain ionospheric sources, such as jet currents. 
There are also other effects produced by the induced 
currents such as changes of phase and of polarization. 
Our point of view in the present discussion is that the 
induced currents are something of a hindrance to our 
investigation of the primary sources of the phenomena, 
and we study them merely in order to get rid of their 
effects. There is of course another point of view with 
regard to these earth currents, namely, that their 
study may give some information about the distribu- 
tion of electrical conductivity within the earth. Two 
methods have been developed for obtaining such infor- 
mation, (i) the magnetic potential method based on 
analyses of magnetic variations over the earth's surface 
or over some suitable part of it, and (ii) the magneto- 
telluric method based on the relations at different 
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frequencies of mutually orthogonal horizontal com- 
ponents of the electric and magnetic variations at a 
given station. A valuable and extensive discussion 
of the theory of the magnetotelluric method has been 
given by Wait [1962]. In the present paper, however, 
we shall be concerned only with the effects of induced 
earth currents on the magnetic field. 

2. Order of Magnitude of the Conductivities 

Involved 

The conductivity of dry surface rocks is frequently 
of order o-= lO'^'mho/mllO'^'^emu) or less, that of 
moist earth and sedimentary rocks 10~^ — 10~2mho/m 
and of sea water, 4 mho/m. Thus the conductivity 
of layers near the surface may vary by several orders 
of magnitude. 

± ii\^ \.^oiiiiiciica ui iiic *^uiiuu\JLiviLy cii ^icaici UcpiiiS) 

have been obtained mainly from analyses of the slower 
magnetic variations, combined with theoretical studies 
of electromagnetic induction. A summary of results of 
these studies is given in figure 1, taken from a recent 
paper by Eckhardt, Larner, and Madden [1963]. It 
will be seen that the mean conductivity near the sur- 
face is estimated at between 10"^ and 10~^ mho/m 
and rises sharply in the region of about 600 km depth 
to over 1 mho/m. The Cantwell-McDonald curve 
shows a rather rapid rise of o" at a few tens of kilo- 
meters depth to a value greater than 5.10"'^ mho/m. 
This part of the profile was deduced by Cantwell from 
magnetotelluric measurements. It should be stated 
that the magnetic variations results obtained by Lahiri 
and Price [1939] certainly did not preclude a high con- 
ductivity at say 60 to 70 km depth, but they indicated 
that the conductivity must in that case decrease 
to appreciably less than 10~*^ mho/m beyond this 
depth and then rise again in the region of about 500 to 
600 km depth. 
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Figure 1. Theoretical electrical conductivity profiles for the 
earth's mantle (after Eckhardt, Larner, and Madden). 



3. Calculation of the Induced Currents 

In the present discussion we consider only the range 
of frequencies from 1 to 0.001 c/s, i.e., periods T from 
1 to 1000 sec. For these values of T (or any greater 
values) and for linear dimensions of the conductor not 
greater than those of the earth, it is permissible to 
ignore displacement currents and use the diffusion 
equation for the electromagnetic field. For values 
of T appreciably less than 1 sec, this is not permissible, 
but the problems which then arise are of a different 
kind, and in these problems it is often sufficient 
to treat the earth as a perfect conductor. 

We shall thus consider only the induction of currents 
by varying magnetic fields. Reference has sometimes 
been made to the possibihty of electrostatic induction 
of currents by varying atmospheric electric fields, 
but apart from thunderstorm effects, any earth 
currents likely to be produced in this way would be 
some orders of magnitude smaller than those produced 
by varying magnetic fields, and will therefore be 
ignored. 

In view of the considerable range of conductivities 
of the crustal layers of the earth, it is clear that those 
induced currents that are near the surface will vary 
greatly in strength from one place to another. Hence 



if the frequency of the primary field is high enough to 
confine the induced currents, because of the "skin 
effect" to the crustal layers, the effects of these 
currents on the observed micropulsation fields will 
differ greatly from one place to another. It is there- 
fore necessary to examine the depth of penetration of 
the induced currents. Before doing this we may note 
that the mathematical problems which arise in the 
study of earth currents are of two kinds. In the first 
we are concerned with inducing fields of global dimen- 
sions and with averaged or smoothed conductivities 
of the earth as a whole; these we may call global prob- 
lems. In the second kind the inducing field may or 
may not be of global dimensions but we are now con- 
cerned with strictly local values of the conductivity, 
we may call these local problems. The problems 
arising in the interpretation of surface observations of 
rapid fluctuations such as micropulsations are often 
of the second kind, but there are exceptions. For 
example, studies of the overall effect of the oceans lead 
to problems of a more global character. On the other 
hand, the investigation of coastal effects is a local 
problem. Local problems are concerned with the 
great differences of conductivity of different geological 
strata. In these local problems we can often ignore 
the sphericity of the earth and treat it as a semi- 
infinite body having a nonuniform distribution of 
conductivity, our attention being confined to some 
limited region. On the other hand, the inducing field 
will often be of large dimensions and therefore effec- 
tively uniform over the region being studied. This 
has sometimes been taken to imply that the actual 
finite dimensions and distribution of the inducing 
field can be ignored in studying the induced currents, 
but this is not always the case. A knowledge of the 
inducing field in the immediate neighborhood of some 
particular station and of the surrounding ground con- 
ductivity is not usually sufficient to determine the 
strength of the induced currents in that neighborhood. 
This will be affected by the nature and distribution of 
the entire inducing field and by the average properties 
of the conductor over a region of corresponding dimen- 
sions. It is also true, of course, that the entire in- 
duced current system will contribute to the magnetic 
field at a particular station, but this fact is not so im- 
portant because usually the largest contribution comes 
from currents in the immediate neighborhood. The 
really important point is that, not only the local 
properties, but the properties of the conductor as a 
whole — its geometry and the distribution of its con- 
ductivity—determine the strength of the currents in 
any neighborhood being studied. 

Thus the problem of calculating the contribution of 
the induced field to the observed total field will 
frequently reduce to finding out how, in a particular 
neighborhood, some extended system of induced 
currents is redistributed by local differences of con- 
ductivity. The extended current system will usually 
need to be determined or estimated from results ob- 
tained in some global type of problem, where the earth 
is treated as a sphere with a suitable smoothed distri- 
bution of bulk conductivity. Many problems of this 
type are discussed in the literature. When the 
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extended average current system has been found in 
this way, the determination of the local redistribution 
of these currents is not purely a problem in the dis- 
turbance of steady current flow, because the currents 
are alternating and a skin effect will usually be in 
evidence. 

4. Separation of the External and Internal 
Sources by Analysis of the Surface Field 

It would not be necessary, of course, to actually 
calculate the induced currents for the purpose of 
removing their effects from the observed surface field 
if one could separate the external and internal (in- 
duced current) sources at any instant by analyzing 
directly the observed surface distribution of the field. 
Theoretically, it would be possible to do this by using 
well-known results in potential theory, but in practice 
we do not have a sufficiently detailed knowledge of the 
surface distribution to do it satisfactorily. Hence it 
is necessary to supplement any studies of the surface 
field by deductions drawn from suitably illustrative 
calculations of the induced currents. 

One of the difficulties inherent in the analysis of the 
surface distribution of the field of fairly rapid geomag- 
netic oscillations is that there are two patterns present 
in this distribution. One is usually a small scale 
pattern determined by the distribution of conductivity 
in the surface layers. This conductivity may change 
greatly in a distance of a few tens of kilometers. 
The other is generally a large scale pattern determined 
by the inducing field, which is often of global dimen- 
sions. If one seeks to determine the general dis- 
tribution and overall dimensions of the field from 
observations at stations over a sufficiently extended 
area, it is obviously necessary to have regard to the 
fact that the field at any one station may be consider- 
ably affected by local geological (possibly subterranean) 
conditions. Due allowance must be made for this by 
estimating the effects from local surveys, or otherwise. 
The local irregularities can of course be reduced to 
a minimum by siting the stations in an area as uniform 
as possible. 

A program of simultaneous observations of various 
types of micropulsations over an extended area would 
certainly provide information of the greatest value 
about the fields involved. It is understood that a 
program of this kind is being organized in the United 
States. One of the most elementary pieces of infor- 
mation we require in order to assess the importance of 
induction effects is a knowledge of the overall extent 
of the inducing fields, because this is an important 
factor in determining the intensity and depth of pene- 
tration of the induced currents. 

Other measurements of great value in the study of 
induction effects are those that have been made on 
ice islands in the Arctic Ocean, such as those reported 
by Hessler [1962 see also Swift and Hessler, 1964], 
Heirtzler [1963], and Zhigalov [I960]. The ocean 
affords a horizontally uniform environment, with the 
conductivity of the sea water known and that of the 
seabed probably very much smaller. Also an ice 
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Figure 2. Curves of the ratio rji^ compared with the ocean 
profile for the ice-island station North Pole 6, February-March 
1958 (after Zhigalov). 

Lat 79°30'-80°44' N, long. 153°57' - 150°17' E. 

1. Ocean depth in meter, 

2. ratio rzlrn, mean for all hours of day, 

3. ratio rzlrn, for variations of period f^ 10 min. 



island, unlike an ordinary oceanic island, does not 
seriously disturb the induced currents in the immediate 
vicinity. Hence the conditions are most favorable 
for interpreting the observations and studying the 
induction effects. Some of the results obtained by 
Zhigalov, from measurements on a drifting ice station 
"North Pole 6" are shown in figure 2, which illustrates 
the attenuation of the vertical component by electric 
currents induced in the sea. The ratios of the ampli- 
tudes of the vertical variations rz to those of the hori- 
zontal variations r// are compared with the depth of 
the sea at a number of positions of the drifting station. 
It will be seen that the ratio rzlrH decreases consider- 
ably as the depth of sea increases, and the curve of 
rzlrH is very similar to the ocean depth profile. Also 
the ratio is considerably smaller for the variations of 
shorterperiod. These results are in accordance with 
the effects expected from currents induced in the sea. 



5. Depth of Penetration of Induced Currents 
and Magnetic Field 

The well-known formula for the skin effect for alter- 
nating currents in a conductor is derived from the solu- 
tion of the diffusion equation 






= cr/jii(oI 



(1) 
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for the current density / in a semi-infinite uniform 
conductor occupying the half-space z < say. The 
period is 27r/a) in seconds, and cr, ijl, are the conduc- 
tivity and magnetic permeability, respectively, in 
MKS units. The permeability /jl may be taken as 
that of free space (fio) in all the cases that need be 
considered. The expression for the current density 
at depth z is then 

/ = /oexp {- zifMoO-iayyl^} (2) 

so that the amphtude contains the attenuation factor 
exp \~zi — ^ — ) [• The skin depth is correspond- 
ingly taken as 



d^O^] 



-1/2 



(3) 



Typical values of d in kilometers, for various conduc- 
tivities and periods are shown in table 1. 



Table 1. Values of the skin depth d in kilometers corresponding 
to various values of the conductivity and period T — 2ttIu) 



Period 


T 


1 


10 


100 


1000 


Conductivity 
in mho/m 


10-^ 

10-2 

4 


50 

5 

0.25 


159 
15.9 
0.8 


500 
50 
2.8 


l.,590 
159 
8.0 



This table suggests that the depth of penetration of 
currents in continental areas is quite considerable 
even for the short periods considered. But it must be 
remembered that the formula used to obtain these 
results is based on a number of simplifying assump- 
tions that may not be apphcable to the geophysical 
problem under consideration. The sphericity of the 
earth is ignored and the conductor is assumed to 
extend downwards to infinity. Also the current flow 
is assumed to be uniform over any horizontal plane. 
It is probable that the assumption that the sphericity 
of the earth can be ignored does not lead to any serious 
error in most cases, but the other two assumptions can 
both conceivably do so, and it is necessary to examine 
them further. 

When the actual distribution and extent of the in- 
ducing field is taken into account, the induced current 
distribution is no longer uniform over each horizontal 
plane, and the depth of penetration of the currents is 
reduced. An estimate of the effect produced by the 
finite scale of the inducing field can be obtained by 
considering the field whose potential, just above the 
surface of the conductor, is 



where P satisfies the equation 



d^P d^P 



(5) 



because H satisfies Laplace's equation. For example, 
if the field is independent of y, P will be a trigonometric 
function, say sin vx, so that the field is periodic in the 
AC-direction with wavelength 27rlv. In the shghtly 
more general case corresponding to (5) we may regard 
27t/v as a length representing the scale of the field 
pattern. 

It can be shown [Price, 1950, 1962] that the density 
at depth z of the resulting induced currents in the con- 
ductor, assumed to extend downwards to infinity, is 
given by 



(6) 



IcTicOlX, JdP dP 

e + v '"" \dy' dx 


where 


d^ = yioi(TO} -\- v'^ . 


The attenuation factor is therefore 


exp \ —z 


71 1 \i/2 1 , 



i;2i 



(7) 



(8) 



a-=-^e^^+^^^PU, y) 



(4) 



which shows that the attenuation is greater than when 
the field is uniform, corresponding to v^O. The 
increase in the attenuation wi ll be si gnificant if v is 
not too small compared with v(/Zo^xo) . 

Now IttIv is not likely to be much less than 4 X lO^m, 
i.e., about four times the height of the E layer where 
some of the current sources may be situated, and it 
cannot be greater than 4 X lO^'m, the circumference of 
the earth. Of course fields of global dimensions can 
only be properly discussed by treating the earth as a 
sphere, but since the depth of penetration is usually 
a small fraction of the earth's radius, it is possible to 
treat the earth's surface as a plane for getting an 
approximate idea of the magnitude of the effect we 
are considering. From the above we see that the 
values of v for natural fields are likely to range from 
1.6 X 10-7 for global fields to 1.6 X lO"'^ for local fields. 

We consider first the case where the conductivity 
is low and the period fairly long so that we have the 
maximum penetration of the field. We take c r= 10" ^ 
mho/m and T= 100 sec. We then have \/(/xoO-w) 
~ 2.8 X 10"^, which is certainly within the likely range 
of values of v. We can infer that, in this case, the 
scale of the inducing field affects the depth of penetra- 
tion appreciably. 

If we take a higher conductivity say o"=10-* 
mho/m — which is probably greater than that of most 
continental strata but still considerably less than that 
of sea water — a nd a s horter period, say F= 10 sec, 
the value of V(/LtoO'a)) becomes 2.8X10"^. Hence, 
in this case, the depth of penetration is not likely to 
be appreciably affected by the scale of any natural 
inducing field. 



1164 



In the case of sea water a ~ 4mho/m, so that 
V(/ioo-a)) ~ 5.7 X 10-3(r)-i/2. Hence, if one could 
treat the sea as of infinite depth, the depth of penetra- 
tion of the induced currents would apparently not be 
sensibly affected by the scale of the inducing field, 
unless the period of the oscillations was at least a few 
hours. But actually the finite depth of the sea has a 
very important effect, and we find that the scale of the 
inducing field becomes important when the depths of 
the real oceans are considered. 



6. Depth of Penetration of Currents and 

Field in the Oceans 

If we consider an ocean of depth D m and assume 
that the conductivity of the seabed is neghgible com- 
pared with that of the sea water itself, the intensity 
of the induced currents at depth z is found to be [Price 
1962] 



M 
m 



o-ez 



{1+/36- 



26(0- z) 



} 



l-|-j3e 



-200 



where 



P = 



e-\-p 



(9) 



(10) 



and is defined by (7). 

The vertical component Z of the field attenuates in 
exactly the same way as the current density /, i.e., 
Z(z)/Z(0) is given by the same formula (9). The hori- 
zontal component H attenuates more rapidly, however, 
and is given by 



H{z) ^ e-^'{1-Pe- 
H{0) l-pe 



2d(D-z) 



-200 



(11) 



These results were obtained from the direct solution 
of the diffusion equation with appropriate boundary 
conditions at z = and z — D. It will be seen that if 
the expressions (9) and (11) are expanded by the 
binomial theorem, they can also be regarded as repre- 
senting the effect at depth z of the incident (evanescent) 
wave e~^^ together with the effects of the sequence of 
waves reflected, with the appropriate transformations, 
from the lower and upper surfaces. The important 
difference between the formulas (9) and (11) should be 
noted. This disappears when D^^o, i.e., all the 
field components attenuate in the same way in a con- 
ductor of infinite depth, but when the depth is finite 
H attenuates more rapidly than / or Z. 

The formulas wiU apply equally well, of course, to 
any highly conducting stratum over a poorly conducting 
basement. More elaborate formulas applicable to 
a series of layers of different conductivities have been 
obtained for purposes of electromagnetic deep sound- 
ing, but these will not be considered here. 

To illustrate the effect of the finite value of D on the 
depth of penetration of the horizontal component of the 
field, the values of the amplitude ratio mod {H{z)IH{0)} 
at various depths z in an ocean of depth 1 km are com- 



pared with the values at the same depths m an ocean 
of infinite depth. The period T is taken as 100 sec and 
the value of v as 10"^. The results are shown in table 2. 



TaB[.E 2. The reduction in amplitude of H at various depths z 
in (i) an ocean of infinite depth, (ii) an ocean of depth 1 km, 
(Hi) an ocean of depth z 









7= 100,.^= 10-* 








z(m) = 


10 


10=^ 


5X10^ 


9X10^ 


103 


2X 103 


5X103 


Z) = oo 


1.0 

0.99 

.30 


0.96 
.90 
.03 


0.82 
.50 
.006 


0.70 
.80 
.004 


0.67 
.003 
.003 


0.45 


0.14 


D = z 


.002 


.001 



In this table the first row of values gives the ratios of 
the amplitude of H at depth z to that at depth zero in 
an ocean of infinite depth, the second row the values in 
an ocean of 1 km depth. The third row gives the values 
at the bottom of the sea when the sea has depth z. 

The remarkable differences in these values should 
be noted. It is clear that the currents are much more 
concentrated near the surface when D is small, and 
the field is therefore much more attenuated than the 
usual skin formula indicates. 

It is also interesting to find that when D is small 
compared with the skin depth, the dimensions of the 
inducing field 27r/^^ become important even when a is 
large. We find, in fact, that when 



D<^ \e\-'<27rlp 



(12) 



the above formula (11) for H{z)IH(0), when z = D, 
reduces to 



1 



H{0) l-\-fjioi(TO)Dlv 



(13) 



which is the well-known formula for shielding by a 
thin conducting plate. Thus doubhng the scale 

— of the inducing field would reduce the amphtude 

V 

ratio for H, i.e., increase the attenuation, bv the same 
amount as doubhng the frequency (w) or the conduc- 
tivity a. 

The above results all relate to the horizontal com- 
ponent H. When the vertical component is con- 
sidered, we find a very different result. Thus for the 
same* values of T, v, as above, the reduction of the 
amplitude of Z at various depths z are shown in table 3. 



Table 3. The reduction in amplitude of Z at various depths z in 
(i) an ocean of infinite depth, (ii) an ocean of depth I km, (Hi) 
an ocean of depth z 









7 = 100,1^=10- 








z(m) 


10 


10* 


5X10* 


9X10* 


103 


2 X IQs 


5X103 


D = oo 


1.00 
1.00 
1.00 


0.% 
.99 
.98 


0.82 
.96 
.93 


0.70 
.91 
.90 


0.67 
.89 
.89 


0.45 


0.14 


D = z 


.82 


.28 
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This table shows that even at the bottom of an ocean 
of depth 2,000 m there is only shght attenuation of 
the amphtude of the vertical component. It must 
however be remembered that the vertical component 
is usually already greatly reduced at the surface by 
the induced currents whereas the horizontal com- 
ponent is increased there. 



7. Effect of Induced Currents on Surface 
Values of H and Z 

The induced currents will tend to reduce the ampli- 
tude of the Z-component and increase that of the 
//-component. For very high frequencies — specifi- 
cally those for which (/xococr)"^ can be taken as zero 
in the calculations — the Z-component will be prac- 
tically extinguished, while, on a flat surfaced conduc- 
tor, H will be doubled. When the curvature of the 
earth has to be taken into account, e.g., when fields and 
conductivity distribution are of global dimensions, 
H would be increased in the ratio (2n+ l)/(7i+ 1), 
where n is the order of the spherical harmonic involved. 

To obtain numerical estimates of the effects for 
other frequencies, we first make calculations for a 
single layer conductor, since, as we have seen, this 
seems an appropriate model for discussing the effects 
of the oceans. 

Denoting the contributions to the surface field 
{H, Z) of the inducing and induced fields by (//o, Zo) and 
(//i, Zi) respectively, we find 



Table 4. The effects of induced currents on H and Z Variations 
at the surface of an ocean of depth D 



Ho Zo l-j82e-2»o' 



(14) 



This formula has been used to estimate the en- 
hancement of the H variations and the reduction of 
the Z variations on the surface of an ocean of depth 
D for various values of D. The scale length 27tIv of 
the inducing field has been taken as lO^m (i.e., 1000 
km), the conductivity o" as 4 mho/m and the period T 
as 1000 sec. The results wiU be unaltered if c and T 
are both changed in the same ratio, e.g., if we take cr 
= 3 mho/m which is likely to be a more representative 
value for the Arctic Ocean, then the results will apply 
for r=750 sec. The results are shown in the first 
two rows of table 4. They show that, for the varia- 
tions considered, the Z component would be practi- 
cally extinguished over an ocean of depth 5 km and 
the H component would be nearly doubled. 

The differences of phase between the inducing and 
induced field, determined from (14) are shown in the 
third row of the table. It will be seen that the phase 
differences are quite small. 

To obtain the relationship between the total Z and 
H variations at the surface we easily find from (14) that 



Z^ Zo + Zi ^Zoi^ l+)8e-^^^ 





^(m) 


5X10^ 


10^ 


2 X lO-' 


3X10^ 


4X103 


5X103 


(1) 


//1///0I 


1.56 


1.71 


1.83 


1.88 


1.91 


1.93 


(2) 


Z./Zol 


0.44 


0.29 


0.17 


0.12 


0.09 


0.07 


(3) 


siTgiH^IH,) 


.0° 


.0° 


. 1° 


.6° 


.9° 


1.1° 


(4) 


rzlrn (calculated) 


.29 


.17 


.09 


.06 


.05 


0.04 


(i3) 


rzjrn (observed) 


.53 


.39 


.15 








(6) 


rzlrn (calculated) 


.37 


.19 


. 10 


.06 






(V) 


rzlrn (calculated) 


.53 


.28 


.14 


.09 







The values in row (4) are calculated from (15) with 27r/i^= 10«, those in row (6) from (16) 
with 2ttIv-=W and those in row (7) from (16) with 2ttIv=1x 101 For all cases o- = 4 
mho/m with T^ 10" sec or, alternatively, o- = 3 mho/m with T= 750 sec. 



The amplitude ratio | ZqJHq \ for the external source 
fields will of course depend on the position of the sta- 
tion where the particular variations are observed. If, 
for simplicity, we assume that the average amplitude 
ratio for these variations is unity, the average ampli- 
tude ratio, say rzlrn^ for the total variations may then 
be calculated from (15). The calculated average val- 
ues of rzlrn are shown in the fourth row of table 4, and 
the values of these ratios estimated from Zhigalov's 
observations (fig. 2) are shown in the fifth row. It will 
be seen that the calculated values are less than those 
obtained from the observations. This could of course 
be due to several different factors. For example, our 
assumption that for the inducing field (rzlrHJo ~ 1, may 
not be true for the variations observed by Zhigalov. 
Again we have taken L to be 1000 km, and the average 
scale factor of the real inducing fields may be greater 
or less than this. Finally, the fact that the ocean has 
a finite area has been ignored, and there is probably 
an important coast effect extending over the shallower 
parts near the edge. 

A rough idea of the nature of the dependence of 
rz/rn on the parameters v, co, o", and D can be obtained 
by treating the ocean as a thin conducting plate of in- 
tegrated conductivity crD. This gives 



Z Zo 



1 



H Ho l-\-iuioi(JO)Dlp 



(16) 



which is, in fact, the approximate form of (15) when 



vD«l 



id \eD\«i. 



(17) 



(15) 



This approximate formula has been used to calculate 
rz/rH for L = 1000 km and L = 700 km, and the results 
are shown in rows 6 and 7 of table 4. It will be seen 
that the latter value of L gives a good fit to the values 
deduced by Zhigalov from the observed amphtude of 
the variations of period ^ 10 min. Though the for- 
mula is not very accurate for the values of \6D\ in- 
volved, we may reasonably conclude that Zhigalov's 
results can be adequately explained by induction of 
currents in the sea. Further, since the induction 
effects appear to be quite sensitive to the scale length 
L (or 27t/v) of the inducing field, it may be possible, 
from studies of similar and more extensive ice island 
observations, to obtain useful estimates of L as well as 
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of the other parameters involved in these induction 
effects. In more elaborate studies, it would, of course, 
be desirable to take into account the possible effects 
of the conductivity and topography of the ocean bottom. 



8. Induction Effects Over Land Areas 

The general agreement between the calculated and 
observed induction effects over the oceans, where 
conditions are fairly uniform, gives promise that, by 
using other and usually rather more elaborate mathe- 
matical models, some idea of the nature and magnitude 
of induction effects over land areas may be obtained. 
A simple model may first be used to get a general idea 
of the magnitude of the effects to be expected. More 
elaborate models must then be used to discuss the 
distribution of these effects and their variation with 
frequency and other parameters. We consider fii^st 
a uniform conductor occupying the half-space z > 0, 
and take various values of the conductivity cr that 
would represent possible average or effective values of 
CT throughout a depth somewhat greater than the skin 
depth for any particular frequency co/27r of the oscilla- 
tions. For this model we find 



we may note that 



Q- ^ ap^TI27TiJLo. 



(19) 



Ho Z() v(/xoicrct) -\-v^)-hp 



(18) 



For the induction effects to be important, /noaw 
must not be small compared with v^, but if it is large 
compared with i^^, then the effects are simply the ex- 
tinquishing of Z and the approximate doubling of H. 
More interesting changes of the induction effects with 
frequency of oscillations will occur when fjLocroj is 
comparable in magnitude with p^. 

Table 5. Values of mod j8 and arg ft for different values of the 
ratio a = iJii)G-(x)lv^ 



a = 


0.01 


0.25 


1.0 


4.0 


25 


100 


inod)Q = 
arg)8 = 


0. 0025 
90° 


0.065 

77° 


0.22 
66° 


0.48 
39° 


0.75 
16° 


0.87 
8° 



Table 5 gives the values of mod j8 and arg y8 for 
values of the ratio a = /jloCtcjIv^. This table shows that 
the induced currents just begin to be effective when a 
reaches 0.25, H then being increased by a factor 1.065 
and Z decreased by a factor 0.935. The effects may 
be regarded as becoming really significant when a 
reaches 1, since H is then increased by 20 percent 
and Z decreased by that amount, so that the corre- 
sponding ratio of the amphtudes of H and Z is then 
two-thirds that for the inducing field alone. When a 
reaches 100, H is increased by a factor 1.87 and Z 
decreased by a factor 0.13. Hence, for a ^ 100, the 
induction effects are predominant and Z is practically 
extinquished compared with H. 

To determine the possible ranges of the parameters 
cr, T, V which correspond to the above values of a 



Now we have already seen that v for the inducing field 
is probably within the range 1.6X10"^ to 1.6X10--^ 
so that 1/2 ranges from 2.5 X lO'^^ for global fields to 
2.5 X 10-1^ for local fields. At the lower hmit of v 
(global fields) with a=l (just significant induction 
effects) we have o" ~ 3 X 10-^ T. This shows that 
induction effects will be significant for these world- 
wide fields for all variations of period r< 10^ sec, 
provided the effective bulk conductivity is greater than 
3 X 10-6 mho/m (3 X lO-^^ emu) and this is probably 
true in almost all cases. It must however be pointed 
out that this result refers only to fields of the widest 
possible extent, i.e., those for which 27rlv = ^X 10^ m, 
and the coefficient a — the "induction response" 
coefficient — varies as v~'^. If we take j^— 1.6X 10-^, 
corresponding to fields extending to about 4,000 km, 
we find from (19) that induction effects for variations 
of period 10^ sec are significant if cr > 3 X lO""* mho/m. 
This value is probably near the upper hmit of the 
effective bulk conductivity. 

The induction effects will be predominant (a =100) 
in fields of maximum global extent for all variations 
whh T< 10*^ sec, if o- > 3 X 10"^ mho/m. For fields 
extending to about 4000 km, and 7=10^ sec, cr would 
have to be greater than 3 X 10"^ mho/m, which seems 
unlikely in most continental areas. Assuming cr is 
of order 3 X 10"^ mho/m, induction effects will pre- 
dominate for these fields when 7< 10 sec. 

For strictly local fields, corresponding io v - 1.6 
X 10-\ cr -^ 3 X 10-^ a T. Hence for T= 10*^ sec, a 
would have to be as high as 3 X 10"^ mho/m for there 
to be significant induction effects, and as high as 3 
mho/m for these effects to be predominant. For 
T= 1 sec, the induction effects are likely to be signi- 
ficant even for these local fields, but will not be pre- 
dominant unless the effective a is greater than 3 X lO"'^ 
mho/m. 

We conclude that for 7 ~ 1 sec, the induction effects 
are always important and will be predominant except 
for inducing fields of very hmited extent. When 
r~10^ sec, the importance of the induction effects 
in continental areas will depend very much on the 
dimensions of the inducing field. 



9. Local Effects Due to Nonuniform 
Conductivity 

The discussion in sections 5-8 has been based on 
the assumption that the earth can be treated as a con- 
ductor with horizontal layers of uniform conductivity. 
In many land areas, however, there are considerable 
variations of cr in the horizontal directions, and in some 
cases actual discontinuities, depending on the geology 
of the area. This nonuniform distribution of con- 
ductivity will affect the distribution of the induced 
currents. Hence the induction effects, particularly 
for short period oscillations, will vary considerably 



1167 



over such an area. The effect of a nonuniform distri- 
bution of conductivity is also shown by the increase of 
rzlrn at land stations near the coastline of a deep ocean 
[Schmucker, 1964]. This increase is caused by the 
induced currents in the relatively highly conducting 
sea water, which tend to concentrate near the boundary 
of the ocean because of self induction. 

The nature of the mathematical problems that arise 
in considering these local features has already been 
referred to in section 3. They are not usually direct 
induction problems of the type considered in the pre- 
vious sections. What is now usually required is a 
method for calculating the disturbance produced by a 
local nonuniform distribution of conductivity in some 
average system — possibly quite extensive — of induced 
currents. In this sense it is a disturbed skin effect 
problem. Very few complete solutions of this type of 
problem, apphcable to oscillations of period 1000 
sec or less, have as yet been published, though a 
number of people are working in this field, and one 
may hope that eventually a hbrary of solutions of 
relevant mathematical problems will be available for 
helping to interpret the observations. 

(Paper 69D8-547) 
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